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\0 ABSTRACT 
O 
O 

^ I A statistical model of the polarization of pulsar radio emission is used to derive 

^ I the general statistics of a polarization vector's orientation angles. The theoreti- 
cal distributions are compared with orientation angle histograms computed from 

I single-pulse, polarization observations of PSR B2020+28. The favorable agree- 
ment between the theoretical and measured distributions lends support to the 

^ , underlying assumptions of the statistical model, and demonstrates, like recent 

\0 \ work on other pulsars, that the handedness of circular polarization is associated 

^ I with the radiation's orthogonally polarized modes. Comprehensive directional 

I statistics of the vector's orientation angles are also derived, and are shown to 

\^ ' follow the Watson bipolar and Fisher distributions in its limiting forms. 

o ■ 
o 



Subject headings: methods: data analysis, statistical - polarization - pulsars: 
general - pulsars: individual (PSR B2020+28) 



1. INTRODUCTION 

X 

' Apart from the classic sweep of the polarization position angle across a pulsar's average 

profile, the orthogonal polarization modes (0PM) are perhaps the most striking feature 
in the polarization of pulsar radio emission. The modes appear in bimodal histograms of 
polarization position angle as two concentrations of data points separated by about 90° (e.g. 
Stinebring et al. 1984). The modes are elliptically as well as orthogonally polarized, so that 
the polarization states of the modes reside at antipodal points on the Poincare sphere (Cordes 
et al. 1978). 0PM are thought to arise from wave propagation effects in the magnetized 
plasma above the pulsar's polar cap (e.g. Allen & Melrose 1982; Barnard & Arons 1986). 
The simultaneous interaction of the modes may be responsible for the depolarization of the 
emission at high radio frequency (Manchester et al. 1975; Morris et al. 1981). 
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The notion that 0PM are eUiptically polarized suggests that the emission's hnear and 
circular polarization should be analyzed together, instead of individually, for a proper inter- 
pretation of pulsar polarization measurements (McKinnon 2003; hereafter M03). Compre- 
hensive methods have been developed recently for displaying and analysing the polarization 
data in this way. For example, Karastergiou et al. (2003) plotted the measured orientation 
angles of a pulsar's polarization vector in a Hammer- Aitoff projection, confirming the associ- 
ation of circular polarization handedness with 0PM in PSR B1133-I-16. Similarly, Edwards 
& Stappers (2004; hereafter ES04) displayed the vector's orientation angles in a Lambert 
equal-area projection to document and explore the intriguing deviations from mode orthog- 
onahty in PSR B0329-I-54. They also showed that the circular polarization in the pulsar is 
generally consistent with an origin in elliptically polarized 0PM. McKinnon & Stincbring 
(1998, 2000) developed a statistical model of superposed 0PM that predicts that polar- 
ization measurements at a given pulse longitude should form a single, prolate ellipsoid in 
the three-dimensional space defined by the Stokes parameters Q, U, and V. The prediction 
is generally supported by the analyses conducted by McKinnon (2004; hereafter M04) and 
ES04, who independently developed eigenvalue techniques to quantify the dimensions and 
orientations of the eUipsoids. While the eigenvalue analyses and the projections of polariza- 
tion vector angles are generally supportive of the superposed 0PM hypothesis, it has yet 
to be demonstrated that the observed orientation angles arc consistent with those predicted 
by the 0PM statistical model. The statistics of the orientation angles derived in M03 are 
limited in their applicability, primarily because they are restricted to the specific case of 
0PM with equal polarization amplitudes. 

The objectives of this paper are to derive a general distribution for a polarization vector's 
orientation angles by expanding the work of M03 and to compare the theoretical distribution 
with the distributions of observed orientation angles. The derivation is presented in §2, and 
the comparison is made in §3. Conclusions are discussed in §4. 

2. THEORETICAL DISTRIBUTIONS OF ORIENTATION ANGLES 

The joint probability density of a polarization vector's amplitude, colatitude, and longi- 
tude must be determined from the Stokes parameters to calculate the individual distributions 
of the orientation angles. From the statistical model of pulsar polarization (M03), the equa- 
tions for the measured Stokes parameters are 



Q = sin 9o cos (t)o{X- 
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U = sin 9o sin 0o(^: 
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V = cos 0o{Xi - X2) + Xn,v (3) 

The random variables Xi and X2 represent the amphtudes of the mode polarization vectors. 
They have means fii and fi2 and standard deviations ai and o"2- In general, Xi and X2 can be 
correlated random variables, but they are assumed to be independent in what follows. The 
Gaussian random variable X-^ accounts for the instrumental noise. It has a zero mean and a 
standard deviation ctn. The fixed angles 9o and 0o are the colatitude and longitude, respec- 
tively, of the diagonal in the Poincare sphere along which the fluctuations in pulsar-intrinsic 
polarization occur. The angles correspond to twice the ellipticity and twice the position 
angle, respectively, of the polarization vector (i.e. 9o — 2xo and 0o = 2V'o in the nomencla- 
ture of ES04). As discussed in M04, the emission may contain additional components, such 
as randomly polarized radiation (RPR). For the purposes of this paper, these components 
are neglected, or alternatively, assumed to be absorbed in the Xn terms provided that their 
fluctuations are Gaussian with a zero mean and that their Q-U-V covariance matrix is a 
scaled identity matrix. 

The basic shape of a data point cluster formed by the equations for the Stokes parameters 
is a prolate ellipsoid. The orientation of the ellipsoid's major axis is determined by 9o and 0o- 
To simplify the derivation of the joint probability density of a polarization vector's amplitude 
and orientation, it is convenient to specify = so that the ellipsoid is symmetric about 
the z-axis in a conventional Cartesian coordinate system. By virtue of this symmetry, the 
longitude of the polarization vector is uniformly distributed over the interval < < 
2tt, and the longitude and colatitude distributions are statistically independent. Since the 
ellipsoid retains its shape on rotation, measured distributions of colatitude and longitude 
can be compared with the theoretical result by determining the orientation of the measured 
ellipsoid's major axis, rotating the Q-U-V data points that form the ellipsoid via matrix 
multiplication so that the major axis of the rotated ellipsoid is aligned with the z-axis, and 
computing histograms of the rotated values of colatitude and longitude. 

When 60 = and when the mode polarization amplitudes, Xi and X2, are independent, 
Gaussian random variables, the joint probability density of the measured polarization vector's 
amplitude and orientation is 



f{r,e,<p) 



sin^ 



<(27r)V2(l + p2)l/2 



exp 



^(1 + p2 sin^ 6) + 2rii cos 9 - ^J? 
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The parameter ^ = — ^2 the mean value of the pulsar intrinsic polarization, and the 
parameter p = (af + a^Y^"^ /a^ is a measure of the intrinsic polarization fluctuations relative 
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to the instrumental noise. Following the procedure outlined in M03, the joint probability 
density can be used to find the distribution of the polarization vector's amplitude. 



a3,[27r(l + p2)]i/2 
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dx (5) 



Similarly, the vector's colatitude distribution is 
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where y is a function of p, s, and 6* given by 
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The colatitude distribution is characterized by two free parameters; the signal-to-noise ratio 
in polarization, s = fi/cr^, and the mode fluctuation ratio, p. The distribution is shown 
in Figure 1 for different values of s and p. The distribution is distinctly bimodal for large 
values of p, and becomes more sharply peaked at large values of s. Equation 6 generalizes 
the case-specific colatitude distributions for s = and p = derived in M03. 

The conditional density of the polarization vector's colatitude, or the colatitude density 
at a fixed value of the polarization amplitude, Vg, is 



f{OK) = 27r 



firo,9,(b) 



sin^ 



/(ro) w(fi;,7) 
where w{k,j) is a normalization factor given by 



exp[«;(7 -I- cos 9)^], 



(8) 



w 



(k, 7)= / exp{K,x^)dx. 
J 7-1 



(9) 



The colatitude conditional density is characterized by two dimensionless quantities; a con- 
centration parameter given by k — r^p^ /2{1 + p'^)cr]^ and a symmetry parameter given by 
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7 = ji/rop^. Examples of the conditional density are shown in Figure 2. The peaks in 
the distribution narrow as k increases. When 7 = (e.g. when /i = 0), the conditional 
density becomes the Watson bipolar distribution (eqn. 38 of M03), which is always sym- 
metric about 9 = Ti/2. As 7 increases, the conditional density is dominated by the strong 
polarization mode and peaks at small values of 9. The parameters k and 7 are not entirely 
independent. For example, when p -C 1, k becomes small and 7 becomes very large, and 
the conditional density becomes the Fisher distribution with a concentration parameter of 
Kf — 27K ~ fxro/cr^ (see eqn. 37 of M03). Statistical treatments of directional data (e.g. 
Fisher et al. 1987) appear to consider the Watson and Fisher distributions as completely 
unrelated. However, the analysis presented here clearly shows that the two distributions are 
special cases of the more comprehensive conditional density given by equation 8. 

Figure 1 shows that the peaks in the colatitude distribution do not occur at the poles 
of the Poincare sphere. In fact, both the locations and the widths of the peaks depend 
upon the values of s and p. The distribution gives the illusion that the data points in the 
Q-U-V cluster form annuli around the poles. The illusion arises from the computation of the 
colatitude distribution, which does not preserve the density of data points in the cluster (see 
also the discussion of equal-area projections in ES04). The polar concentration of data points 
is retained in the distribution of m = cos 9, where the locations of the 0PM peaks occur at 
9 — 0,7r and only the widths of the peaks are affected by s and p. The distribution of u can 
be computed directly from the distribution of 9 given by equation 6. The distribution is 
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where y is given by 
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The conditional density of u can be derived from equation 8. 
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Fig. 1. — Distributions of a polarization vector's colatitude, 9 (eqn. 6). The top panel shows 
colatitude distributions for a fixed signal-to-noise ratio in polarization, s, with varying mode 
fluctuation ratio, p. The bottom panel shows colatitude distributions for fixed p and varying 
s. 
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Fig. 2. — Conditional density of a polarization vector's colatitude (eqn. 8). The distributions 
are shown for different concentration parameters, k, and symmetry parameters, 7. The top 
panel shows distributions for varying 7 with k fixed. The bottom panel shows distributions 
for varying k with 7 fixed. 
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3. DISTRIBUTION COMPARISON 



To compare the measured distributions of a polarization vector's colatitude and longi- 
tude with the analytical results, the Q-U-V data points comprising the polarization ellipsoid 
must be rotated so that the ellipsoid's major axis is aligned with the z-axis. From Fisher et 
al. (1987), the matrix multiplication operation that performs the rotation is 





'Qo' 




"0" 




Uo 
















" cos do cos 0O cos do sin 0o — sin 9o 

— sin0o COS0O Uo = ^ ■ (13) 

_ sin 9o cos 0o sin 9o sin 0o cos 9o 

Care must be exercised in choosing appropriate estimates of the ellipsoid's orientation 
angles, 9o and 0o, before the rotation is performed. The angles could be estimated from the 
mean values of the Stokes parameters at the pulse longitude of interest from 



arctan(-y^), (14) 



9o — arccos 



{Q) 

{V) 



(15) 



which would suffice for locations within the pulse where the polarization signal-to-noise ratio 
is high, s ^ 0. But at locations where s ~ 0, the mean Stokes parameters will be small, and 
consequently equations 14 and 15 will produce poor estimates of 0o and 6o- Alternatively, 
one could compute the covariances of the Stokes parameters, and the directional cosines of 
the principal eigenvector of the 3x3 covariance matrix could be used to estimate 4>o and 9o- 
The polarization eUipsoid can be highly elongated where 0PM occurs, and the directional 
cosines of the principal eigenvector can be well determined, even where s = 0, because its 
eigenvalue, th = clf(l + P^) (M04), is a strong function of the mode fluctuation ratio, p. 
However, this method produces poor estimates of the orientation angles at locations where 
p ~ 0, even if s 0, because the eigenvalues in this case are equal (rn = T22 = T33 = Cat), 
and thus the eigenvectors are not unique. The best estimate of the ellipsoid's orientation 
angles, and the one used in the analysis presented below, comes from the directional cosines 
of the principal eigenvector determined from the second moments of the Stokes parameters. 
The eigenvalue of the second moment principal eigenvector is Th = cr|^(l + + p^), and the 
directional cosines of the eigenvector will be well determined when either s ~ or p ~ 0. 
When both s and p are small, no polarized signal is present, the polarization ellipsoid is a 
spheroid centered on the origin of Q-U-V space, and there is no need to rotate the spheroid 
(i.e. the distribution of data points is isotropic). 
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Histograms of the measured values of colatitude and longitude were computed at differ- 
ent locations within the pulse of PSR B2020+28 using the 1404 MHz single pulse, polariza- 
tion observations of Stinebring et al. (1984) as follows. From the multiple measurements of 
the Stokes parameters taken at each location, the second moments of the Stokes parameters 
were computed and used to find the directional cosines of the polarization ellipsoid's principal 
eigenvector (major axis). The directional cosines were then used to compute the orientation 
angles of the major axis. The resulting values of Og and (j)o were used in equation 13 to 
rotate the Q-U-V data points that formed the polarization ellipsoid. Finally, the rotated 
Stokes parameters were used to calculate corresponding values of colatitude and longitude, 
which were binned into 50 intervals spread over the appropriate range of angles to form the 
histograms. 

The histograms at three different locations in the pulse are compared with the theoretical 
distributions of the polarization vector's colatitude and longitude in Figure 3. The colatitude 
distributions are shown in the left column of panels in the figure, and their corresponding 
longitude distributions are shown in the right column of panels. The solid horizontal line in 
each of the longitude histogram panels represents the uniform distribution expected from the 
0PM statistical model. The generally good agreement between the theoretical and measured 
longitude distributions, combined with the observation that the ratios of the minor axes of the 
Q-U-V data point clusters are no larger than 1.09 at these locations within the pulse (M04), 
suggests that the polarization ellipsoids are rotationally symmetric about their major axes. 
The observed symmetry of the eUipsoids also suggests that the longitude and colatitude of the 
polarization vector are statistically independent, as predicted by the 0PM statistical model. 
However, the symmetry does not rule out the possibility that (f) and 9 can be correlated or 
dependent upon one another. Indeed, the deviations from orthogonality documented in ES04 
and M04 are evidence that the two angles can be correlated. The smooth, solid curves shown 
in the colatitude histogram panels represent the best fits of the data to equation 6. The best 
fit values of p and s are annotated in each panel. The excellent fits show that equation 6 
is capable of representing a broad range of polarization behavior. Furthermore, these fits 
and the measurements of cluster dimensions by M04 together imply that the clusters aren't 
just any eUipsoids, but the ellipsoids expected from the 0PM statistical model, suggesting 
quantitative, in addition to qualitative, consistency between the model and the data. The 
data used to construct the histograms had significant circular polarization, and the ellipsoidal 
shape of the Q-U-V data point clusters illustrates the association of circular polarization 
handedness with 0PM in PSR 2020+28, similar to what was found in PSR 1133+16 by 
Karastergiou et al. (2003) and in PSR B0329+54 by ES04. Only one polarization mode is 
apparent (p = 0) in the colatitude histogram shown in the top left panel of the figure. In the 
lower panels, the other polarization mode becomes more obvious as indicated by the bimodal 
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histograms and the increasing value of p. Additionally, the values of s in the panels decrease 
as p increases, which is to be expected since the simultaneous interaction of the orthogonal 
modes depolarizes the emission. The fitted values of s are low when compared to what is 
expected from the mean polarization and the off-pulse instrumental noise. However, the 
sizes of the ellipsoids' minor axes, and thus the effective noise of the observations, are much 
larger than the off-pulse noise, causing the lower-than-expected values of s. An emission 
component such as RPR could contribute to the large effective noise (M04, ES04). 

The measured and theoretical distributions of cos 9 arc compared in Figure 4. The data 
shown in each panel of the figure arc the same data used in the corresponding panel of the 
left column of Figure 3. The histogram in each panel is the measured distribution of cos 6'. 
The smooth curve is the theoretical distribution computed from equation 10 using the same 
values of p and s as in Figure 3. In all panels, the data are concentrated at one or both of 
the poles. The agreement between the theoretical and measured distributions is excellent. 

4. DISCUSSION AND CONCLUSIONS 

The statistical model of pulsar polarization developed by McKinnon & Stinebring (1998) 
predicts that multiple measurements of the variable Stokes parameters at a fixed pulse lon- 
gitude will form a single, prolate ellipsoid in Poincare space. The analyses presented in ES04 
and M04 show that Q-U-V data point clusters in PSR B0329+54, PSR B1929+10, and PSR 
B2020-I-28 are indeed prolate ellipsoids. In this paper, the favorable comparison between 
the theoretical and measured distributions of polarization vector colatitude and longitude 
in PSR B2020+28 suggest that the data point clusters are consistent with the ellipsoid ex- 
pected from the model. This implies quantitative consistency between the model and the 
data, and supports the fundamental, underlying assumption of the statistical model which is 
the observed polarization is determined by the simultaneous interaction of two, orthogonally 
polarized modes. 

The analysis of orientation angles presented here and the observations and analyses of 
Cordes et al. (1978), Karastergiou et al. (2003), and ES04 clearly demonstrate the associa- 
tion of circular polarization handedness with 0PM, confirming that the modes are elliptically 
polarized and occupy antipodal points on the Poincare sphere. This result could be inter- 
preted to conclude that the circular polarization is formed by the appropriate phasing of 
the naturally occurring, linearly polarized modes through mode coupling (e.g. Lyubarskii 
& Petrova 1998; Petrova 2001) as the radiation propagates through the pulsar's magneto- 
spheric plasma. If this were the case, we would only observe circular polarization in parts 
of the pulse where both modes were found. However, there are clear examples of circular 
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Fig. 3. — Comparison between the measured and theoretical distributions of a polarization 
vector's colatitude and longitude at three different locations in the pulse of PSR B2020+28. 
The longitude histograms are generally consiscnt with the uniform distributions predicted 
by the statistical 0PM model, as indicated by the solid horizontal line. The smooth curves 
through the colatitude histograms are best fits to equation 6. The fitted values of signal- 
to-noise ratio, s, and mode fluctuation ratio, p, are annotated in each colatitude histogram 
panel. 
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cos(0) 



Fig. 4. — Comparison between the measured and theoretical distributions of cosine colati- 
tude, cos e, at the three locations in PSR B2020+28. 
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polarization where only one mode occurs (e.g. top left panel of Fig. 3 and Fig. 3 of M04). 
So, unless one mode just so happens to be completely "coupled" in the process of creating 
circular polarization, other propagation effects, such as the gyrotropy of the magnetospheric 
plasma, or cyclotron absorption (Melrose 2003) must be responsible for the production of 
circular polarization in the radio emission. 

Karastergiou et al. (2003) and ES04 displayed the orientation angles of the polarization 
vector at a fixed pulse longitude using Hammer-Aitoff and Lambert equal-area projections. 
The colatitude and longitude histograms computed in §3 arc simply another method for con- 
veying the same information. Both the projections and the histograms effectively illustrate 
the fluctuations in the colatitude and longitude of the emission's polarization vector. 

A comprehensive version of the conditional density for the polarization vector's orien- 
tation angles was derived. In its limiting forms, the conditional density follows the Watson 
bipolar distribution when the polarization amplitude is small (/x = 7 = 0) and follows the 
Fisher distribution when the polarization fluctuations are small compared to the effective 
noise (p <^ 1). Since both the Watson and Fisher distributions have many applications in 
the statistical analysis of scientiflc data (M03), the comprehensive version of the conditional 
density given by equation 8 may have more widespread applications. 

The small values of polarization signal-to-noise ratio found in the orientation angle 
analysis and the inflated Q-U-V data point clusters found by M04 and ES04 are reminders 
that the effective noise of single pulse, polarization observations is much larger than the noise 
measured off the pulse. The additional noise may be due to RPR. The colatitude distribution 
(eqn. 6) was derived assuming Gaussian fluctuations in mode polarization amplitudes, and it 
is remarkable that the equation fits the data so well, regardless of the actual type of pulsar- 
intrinsic polarization fiuctuations. Perhaps the good match between equation and data can 
be attributed to the idea that the observed fiuctuations are heavily infiuenced by the large, 
effective Gaussian noise. The caution to interject here is that we are sifting through a large 
effective noise component in analyzing these data, irrespective of how low the off-pulse noise 
might be. Consequently, it is important to account for the effective noise in statistical models 
so that the data can be properly interpreted. 

Again, I am indebted to Dan Stinebring for his continued generosity in providing the 
data used in the analysis. I thank an anonymous referee for comments that improved the 
manuscript. 
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